NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


FLOW  CHARACTERISTICS  OF  A  MULTIPLE  NOZZLE 
EXHAUST  GAS  EDUCTOR  SYSTEM 

by 

Dennis  Leo  Ryan  III 
March  1981 

Thesis  Advisor:  P.F.  Pucci 

Approved  for  public  release;  distribution  unlimited. 


8l  9  11  018 


SKCUMITV  CL AtMriCATION  0*  TNIt  P  A<1»  ftlW  Dm*  «HHW< 


rcrHnrrnrrr  ran  jits 


i 


Item  20:  continued. 

"-'in  terms  of  non-dimensional  parameters  governing  the  flow 
phenomena  from  a  one  dimensional  analysis  of  a  simple 
eductor  system.  The  eductor  system’s  pumping  capacity  was 
improved  over  that  of  the  previously  tested  unmodified 
mixing  stack. 


DD  Form  1473 
V  d  TOjmm-moi 


2 


Unclassifed 

MCUMtV  ChMftMCAVWI  •'  T»*  *»•( 


n 


Approved  for  public  release;  distribution  unlimited 


Flow  Characteristics  of  a  Multiple  Nozzle 
Exhaust  Gas  Eductor  System 


by 


Dennis  Leo  Ryan  III 

Lieutenant  Commander,  United  States  Navy 
B.S.,  United  States  Naval  Academy,  1970 


Submitted  in  partial  fulfillment  for  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
March  1981 


Chairman,  Department  o 


nicalEagineering 


Dean  of  Science  and 


Engineering 


3 


ABSTRACT 


Cold  flow  tests  of  a  four  nozzle  eductor  system  were 
conducted  to  evaluate  the  flow  characteristics  of  several 
mixing  stack  configurations.  A  previously  tested  mixing 
stack  design  used  a  plug  to  shield  the  primary  flow  nozzles 
from  view.  Flow  visualization  was  used  to  determine  the 
flow  pattern  in  the  stack.  The  results  of  the  visualization 
runs  were  then  used  to  modify  the  mixing  stack  geometry. 

The  eductor  system  flow  characteristics  were  evaluated  in 
terms  of  non-dimensional  parameters  governing  the  flow 
phenomena  from  a  one  dimensional  analysis  of  a  simple  eductor 
system.  The  eductor  system's  pumping  capacity  was  improved 
over  that  of  the  previously  tested  unmodified  mixing  stack. 


TABLE  OF  CONTENTS 


I.  INTRODUCTION  .  14 

II.  EDUCTOR  ANALYSIS  .  19 

III.  EXPERIMENTAL  CORRELATION  .  20 

IV.  MODEL  GEOMETRIES  .  21 

A.  MIXING  STACK  WITH  PORTED  AND 

SHROUDED  PLUG . - .  21 

B.  CUSP  MODIFICATION  OF  THE  PLUG  .  22 

C.  MIXING  STACK  SIDE  MODIFICATION  .  23 

V.  EXPERIMENTAL  FACILITY  . 24 

A.  PRIMARY  AIR  SYSTEM .  24 

B.  SECONDARY  AIR  PLENUM  .  26 

C.  TERTIARY  AIR  PLENUM .  2  7 

D.  ALIGNMENT  .  29 

E.  INSTRUMENTATION  - .  29 

VI.  EXPERIMENTAL  METHOD  .  32 

A.  PUMPING  COEFFICIENTS  .  32 

B.  INDUCED  AIR  FLOW .  33 

C.  PRESSURE  DISTRIBUTIONS  .  34 

D.  EXIT  VELOCITY  PROFILES  .  34 

E.  FLOW  VISUALIZATION  .  35 

VII.  DISCUSSION  OF  EXPERIMENTAL  RESULTS  .  37 

VIII.  CONCLUSIONS  .  42 

IX.  RECOMMENDATIONS  FOR  FURTHER  STUDY  .  43 

5 


LIST  OF  REFERENCES  .  44 

FIGURES  .  4  5 

TABLES  . - .  99 

APPENDIX  A:  ONE  DIMENSIONAL  ANALYSIS  OF  A 

SIMPLE  EDUCTOR . - .  140 

APPENDIX  B:  FORMULA  .  152 

APPENDIX  C:  DATA  ACQUISITION  AUTOMATION  .  155 

APPENDIX  D:  UNCERTAINTY  ANALYSIS  .  164 

INITITAL  DISTRIBUTION  LIST  .  166 


6 


LIST  OF  FIGURES 


FIGURE  1  Eductor  Model  Testing  Facility  .  45 

FIGURE  2  Eductor  Model  Testing  Facility,  Secondary 

and  Tertiary  Air  Plenums  .  46 

FIGURE  3  Schematic  of  Mixing  Stack  Eductor  .  47 

FIGURE  4  Dimensional  Illustration  of  Mixing 

Stack  with  Symmetric  Plug  .  48 

FIGURE  5  Schematic  of  Mixing  Stack  With 

Ported  and  Shrouded  Plug  - . .  49 

FIGURE  6  Schematic  of  Plug  Porting 

Arrangement  .  50 

FIGURE  7  Cusp  Illustration  .  51 

FIGURE  8  Dimensional  Illustration  of  Mixing 

Stack  With  Cusp  . - . - .  52 

FIGURE  9  Mixing  Stack  with  Cusp  Installed .  53 

FIGURE  10  Dimensional  Illustration  of  Mixing 

Stack  With  Cusp  and  Side  Modification  .  54 

FIGURE  11  Mixing  Stack  With  Cusp  and  Side 

Modification . .  55 

FIGURE  12  Secondary  and  Tertiary  Air  Plenums  .  56 

FIGURE  13  Tertiary  Air  Plenum  .  57 

FIGURE  14  Mixing  Stack  Entrance  With  Air  Seal  .  58 

FIGURE  15  Mixing  Stack  Exit  With  Air  Seal  .  59 

FIGURE  16  Mixing  Stack  Mounting  Stand  .  60 

FIGURE  17  Mixing  Stack  With  Ported  Shroud  .  61 

FIGURE  18  Schematic  of  Instrumentation  .  62 

FIGURE  19  Schematic  of  Instrumentation  for 

Primary  Air  Flow  Measurement  .  63 


7 


FIGURE 

FIGURE 

FIGURE 

FIGURE 

FIGURES 
a,  b, 

FIGURES 

a-g 

FIGURES 

a-1 

FIGURES 
a-  i 

F?GURE 

FIGURE 

FIGURES 
a ,  b 


1 


20  Velocity  Profile  Apparatus  .  64 

21  Example  Plot  of  Pumping  Coefficient  . 65 

22  Schematic  of  Velocity  Traverse 

Directions  . 66 

23  Schematic  of  Flow  Visualization  Results  ----  67 


24  Experimental  Pumping  Coefficient 

c  Comparisons  . - . 

25  Experimental  Film  Pumping 

Coefficient  Comparisons  . 

26  Axial  Pressure  Distribution 

Comparisons  . 

27  Velocity  Traverses  . 


A-I  Simple  Single  Nozzle  Eductor 

System  . . 151 

C-I  Schematic  of  Data  Acquisition 

System  . 161 

C-II  Illustration  of  Data  Acquisition . 162 


8 


LIST  OF  TABLES 


TABLES  I  through  XVIII 

Computer  Output  of 
Experimental  Data  . 

. 99 

TABLE  C-I 

Sample  Calibration 
Program  . 

. . 160 

9 


NOMENCLATURE 


F 

» 


English  Letter  Symbols 


A 

c 

C 

D 

F 


fr 

’c 


h 

k 

L 

P 


Area,  in.2 

Sonic  velocity,  ft/sec 

Coefficient  of  discharge 

Diameter,  in. 

Thermal  expansion  factor 

Wall  skin- friction  force,  Ibf 

Proportionality  factor  in  Newton's  Second 
Law,  gc  =  32.174  lbm-ft/lbf -sec2 

Enthalpy,  Btu/lbm 

Ratio  of  specific  heats 

Length,  in. 

Pressure,  in.  H  0 


P 

a 

P 

v 

PMS 

R 

s 

S 

T 

u 

U 


2 

Atmospheric  pressure,  in.  Hg 

Velocity  head,  in.  H  0 

2 

Static  pressure  along  length  of  mixing 

stack  in.  H  0 
2 

Gas  constant  for  air,  53.34  ft-lbf/lbm-R 

Entropy,  Btu/lbm-R 

Primary  dimension  of  mixing  stack 

Absolute  temperature,  R 

Internal  energy,  Btu/lbm 

Velocity,  ft/sec 


10 


V 


Specific  volume,  ft3/lbm 
W  -  Mass  flow  rate,  lbm/sec 

D  -  Distance  from  primary  nozzle  exit  to  mixing 

stack,  in. 

Y  -  Expansion  factor 


Dimensionless  Groupings 


AR 

f 

K 

K 

e 

Km 


Secondary  flow  area  to  primary  flow  area 
ratio 

Area  ratio 
Friction  factor 
Flow  coefficient 

Kinetic  energy  correction  factor 

Momentum  correction  factor  at  the  mixing 
stack  exit 

Momentum  correction  factor  at  the  primary 
nozzle  exit 


M 


Mach  number 


AP*  =  p* 
PMS* 

Re 

D/S 

T* 

T£  =  TT* 
W*  =  W* 


Pressure  coefficient 

Mixing  stack  pressure  coefficient 

Reynolds  number 

Standoff;  Ratio  of  distance  from  primary 
nozzles  to  entrance  of  mixing  stack,  (D) 
to  primary  dimension  of  mixing  stack  (S) 

Absolute  temperature  ratio  of  the  secondary 
flow  to  primary  flow 

Absolute  temperature  ratio  of  the  tertiary 
flow  to  primary  flow 

Secondary  mass  flow  rate  to  primary  mass 
flow  rate  ratio 


11 


1 


W*  =  WT* 


Tertiary  mass  flow  rate  to  primary  mass 
flow  rate  ratio 


L/S 

o* 


Ratio  of  distance  from  entrance  of  mixing 
stack  to  diameter  of  mixing  stack 

Induced  flow  density  to  primary  flow  density 


Greek  Letter  Symbols 

U  -  Absolute  viscosity,  lbf-sec/ft2 

p  -  Density,  lbm/ft3 


Subscripts 

0 

1 

2 

f 

m 

or 

P 

s 

t 

u 

w 


Section  within  secondary  air  plenum 

Section  at  primary  nozzle  exit 

Section  at  mixing  stack  exit 

Film  or  wall  cooling 

Mixed  flow  or  mixing  stack 

Orifice 

Primary 

Secondary 

Tertiary  (Cooling) 

Uptake 

Mixing  stack  inside  wall 


Computer  Tabulated  Data  and  Illustrative  Plots 
UPTMACH  -  Uptake  Mach  number 

PA-PS  -  Static  pressure  at  mixing  stack  entrance, 

in.  H  0 
2 

Static  pressure  in  tertiary  air  plenum, 

in.  H  0 
2 


PA-PT 


I.  INTRODUCTION 


With  gas  turbines  becoming  a  more  popular  means  of  power¬ 
ing  naval  vessels,  special  considerations  need  to  be  given  to 
their  particular  air  breathing  and  exhausting  characteristics. 
With  air- fuel  ratios  of  four  to  five  times  that  of  conventional 
steam  plants  and  the  requirement  for  a  relatively  large  amount 
of  combustion  ai",  a  large  quantity  of  hot  exhaust  gas  is 
generated.  Due  to  gas  turbine  design,  these  exhaust  gases 
are  a  temperatures  significantly  above  those  of  conventionally 
powered  ships.  A  few  of  the  problems  caused  by  these  high 
temperatures  are  thermal  damage  to  electronic  equipment  located 
on  the  masts  of  these  ships,  hot  gas  corrosion  of  masts  and 
other  superstructures  located  in  the  hot  gas  wake,  and  a  sig¬ 
nificant  infrared  radiation  signature  created  by  the  hot  gas 
plume  and  hot  external  surfaces  of  the  stack. 

This  thesis  is  an  extension  of  research  done  by  Ellin 
[Ref.  1],  Moss  [Ref.  2],  Lemke  and  Staehli  [Ref.  7],  and 
Shaw  [Ref.  8]  to  determine  better  geometric  designs  for  the 
exhaust  plenum  and  mixing  stack  system  of  gas  turbine  powered 
naval  ships. 

Ellin  initiated  the  work  by  constructing  an  eductor  model 
testing  facility  consisting  of  an  uptake,  primary  flow  nozzle, 
mixing  stack,  a  means  to  control  and  measure  the  primary  air 
flow,  and  a  means  to  measure  the  secondary  air  flow;  see 


14 


*3 


Figures  1  and  2.  The  primary  air  flow  in  the  testing  facility 
represents  a  gas  turbine's  hot  exhaust  gas.  The  secondary  air 
flow  is  ambient  air  induced  into  the  entrance  of  the  mixing 
stack  by  the  primary  air  flow;  see  Figure  3.  From  Ellin's 
study  of  multiple  nozzle  flow  systems  consisting  of  several 
identical  round  nozzles,  it  was  determined  that  four  primary 
flow  nozzles  were  preferable  to  either  three  or  five,  and 
that  nozzle  length  has  little  or  no  effect  on  the  eductor 
system's  overall  performance.  Ellin  then  verified  the  inde¬ 
pendence  of  the  one-dimensional  model  correlation  parameters 
used  on  flow  rate  or  Mach  number.  He  determined  that  for 
Mach  numbers  from  50%  to  145%  of  the  design  Mach  number  of 
0.064,  the  correlation  parameters  suggested  in  the  one-dimen¬ 
sional  analysis  did  in  fact  provide  a  good  correlation  of  the 
data. 

Moss'  work  followed,  and  it  initially  consisted  of 
verifying  the  one-dimensional  analysis  as  did  Ellin.  He  then 
tested  the  effect  of  the  stand-off  distance  (that  distance 
between  the  exit  plane  of  the  primary  flow  nozzles  and  the 
entrance  plane  of  the  mixing  stack).  For  the  primary  flow 
nozzles  he  tested,  Moss  determined  that  the  optimum  stand-off 
distance  for  maximum  eductor  pumping  was  a  distance  equal  to 
0.5  diameters  (0.5Dm)  of  the  mixing  stack.  An  independent 
investigation  of  this,  conducted  by  Harrel  [Ref.  3],  confirmed 
Moss'  findings.  Moss  then  investigated  the  effects  of  a 
conical  transition  placed  on  the  entrance  to  the  mixing  stack. 
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He  concluded  that  a  straight  mixing  stack  without  an  entrance 
transition  provided  a  better  system  performance. 

The  study  conducted  by  Lemke  and  Staehli  [Ref.  7]  inves¬ 
tigated  the  effects  on  the  eductor  system's  overall  performance 
of  varying  the  geometric  configuration  of  the  mixing  stack 
and  changing  the  area  of  the  primary  flow  nozzles.  Their 
work  showed  that  a  decrease  in  the  ratio  of  the  area  of  the 
mixing  stack  to  primary  flow  nozzles  from  3.0  to  2.5  decreased 
uptake  back  pressure  but  reduced  the  pumping  coefficient  of 
the  eductor.  Lemke  and  Staehli  then  investigated  the  effects 
of  a  solid  diffusor,  a  two-ring  and  a  three-ring  diffuser. 

The  results  of  these  tests  showed  a  decrease  in  uptake  back 
pressure  and  a  improvement  in  the  eductor's  pumping  capacity. 
They  then  performed  tests  on  a  ported  mixing  stack.  Their 
work  determined  that  significant  air  flow  through  the  ports 
could  provide  film  cooling  on  the  inside  of  the  mixing  stack. 

To  enhance  the  film  cooling  provided  by  the  ported  mixing 
stack,  Lemke  and  Staehli  placed  a  shroud  around  the  mixing 
stack.  This  shroud  did  not  degrade  the  pumping  or  mixing 
characteristics >  yet  it  provided  thermal  shielding  of  the 
mixing  stack.  Their  final  configuration  was  a  combination  of 
the  ported  mixing  stack  with  flow  through  shroud  and  diffusor. 
Lemke  and  Staehli  concluded  that  this  geometric  configuration, 
combined  with  a  ratio  of  the  area  of  the  mixing  stack  to  the 
primary  flow  nozzles  3.0,  provided  the  best  system  overall 
performance. 


The  work  done  by  Shaw  [Ref.  8]  concentrated  on  a  method 
to  further  reduce  the  infrared  signature  of  the  mixing  stack. 
Of  particular  concern  was  the  concealment  of  the  hot  primary 
flow  nozzles  from  an  overhead  view.  Using  the  results  of 
Lemke  and  Staehli,  Shaw  designed  a  rectangular  mixing  stack 
with  an  internal  symmetric  plug.  Figure  4.  The  internal  plug, 
which  shields  the  primary  nozzles  from  direct  view,  is  cooled 
by  two  means.  The  plug  was  ported,  to  provide  film  cooling 
along  the  exhaust  gas  side  of  the  device.  Also,  the  plug  was 
shrouded  (on  the  non-exhaust  gas  side)  to  direct  the  flow  of 
cooling  air  to  the  ports  and  provide  a  means  of  convectively 
cooling  the  downstream  half  of  the  plug,  see  Figure  5.  Shaw 
found  that  the  rectangular  mixing  stack  with  plug  increased 
the  system  back  pressure  and  reduced  the  pumping  capacity  of 
the  eductor.  He  also  determined  that  the  design  and  place¬ 
ment  of  the  plug  improved  mixing  of  the  exhaust  gases  with 
the  secondary  cooling  air  and  that  the  plug  provided  a 
potential  means  of  presenting  a  cool  surface  to  view. 

Based  upon  the  conclusions  and  recommendations  of  Shaw, 
it  was  decided  that  geometric  changes  in  the  eductor  should 
be  investigated  in  order  to  lessen  the  penalties  involved  with 
plug  installation.  These  changes  will  be  described  in  detail 
in  Section  IV. 

Evaluation  of  eductor  system  performance  was  measured  in 
four  areas:  the  amount  of  secondary  air  flow  induced  by  the 
primary  air  flow,  the  degree  of  mixing  of  primary  and  induced 
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air  flows  within  the  mixing  stack  system,  the  amount  of  uptake 
back  pressure  impressed  upon  the  turbine  exhaust  by  the  eductor 
system,  and  the  amount  of  wall  cooling  air  available  to  reduce 
the  exterior  stack  temperature  of  the  eductor  system.  Addi¬ 
tionally,  a  qualitative  indication  of  the  flow  patterns  in  the 
mixing  stack  was  obtained  through  flow  visualization. 

The  key  factor  which  allows  cold  flow  testing  to  predict 
the  effects  of  a  hot  gas  eductor  system  is  the  similarity  of 
the  momentum  and  energy  transfer  mechanisms  in  turbulent  flows. 
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II.  EDUCTOR  ANALYSIS 


This  investigation,  being  an  extension  of  the  work  of 
Shaw  [Ref.  8],  and  others  [Ref.  1,  2,  and  7];  uses  the  same 
one-dimensional  analysis  to  model  the  eductor  system.  This 
one -dimensional  analysis  is  given  in  Appendix  A.  As  a  result 
of  this  analysis,  dimensionless  parameters  controlling  the 
phenomena  are  obtained  which  are  then  used  in  the  correlation 
of  experimental  data.  Three  such  dimensionless  groups  used 


a  pressure  coefficient  which 
compares  the  pumped  head  Pa  -PQs 

to  the  driving  head  U  *  of  the 
primary  flow  p 

a  flow  rate  ratio,  secondary  to 
primary  mass  flow  rate 


an  absolute  temperature  ratio, 
secondary  to  primary 
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EXPERIMENTAL  CORRELATION 


For  the  geometries  and  flow  rates  investigated,  it  was 
confirmed  by  Ellin  and  Moss  that  a  satisfactory  correlation 
of  the  variables  P*,  T*,  and  W*  takes  the  form 

_AP1  =  0 (W*T*n)  (1) 

where  the  exponent  n  is  determined  to  be  equal  to  0.44.  The 
details  of  the  determination  of  0.44  as  the  correlating  ex¬ 
ponent  for  the  geometric  parameters  of  the  models  tested  is 
given  in  Reference  [1].  To  obtain  an  eductor  model's  pumping 
characteristic  curve,  the  experimental  data  is  correlated  and 
analyzed  using  equation  (1)  that  is  P*/T*  is  plotted  as  a 
function  of  This  correlation  is  used  to  predict 

the  open  to  the  environment  operating  point.  Variations  in 
the  eductor  model's  geometry  will  change  the  appearance  of 
the  pumping  ability  between  models.  For  ease  of  discussion, 
W*T*°’^4  will  henceforth  be  referred  to  as  the  pumping  co¬ 
efficient.  Similarly,  WT*TT*^'^  will  be  referred  to  as  the 
film  cooling  or  tertiary  pumping  coefficient. 
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IV.  MODEL  GEOMETRIES 


The  multiple  nozzle  eductor  systems  studied  here  are  de¬ 
signed  specifically  for  service  onboard  gas  turbine  powered 
ships.  The  model  consisted  of  a  single  primary  uptake,  a 
single  cluster  of  four  primary  nozzles  of  constant  cross 
section,  and  a  mixing  stack. 

A.  MIXING  STACK  WITH  PORTED  AND  SHROUDED  PLUG 

The  eductor  mixing  stack  as  designed  by  Shaw  with  its 
tertiary  air  flow  cooling  ports  and  associated  plug  shroud¬ 
ing  provided  the  base  design  for  the  initiation  of  this 
study.  Shaw's  design  utilized  a  symmetric  plug  to  shield 
the  hot  primary  exhaust  nozzles  from  direct  view.  The  ini¬ 
tial  design  by  Shaw  is  illustrated  in  Figure  4.  He  then 
provided  two  means  of  cooling  the  surfaces  of  the  plug  that 
were  exposed  to  direct  view.  The  first  method  provides  film 
cooling  along  both  sides  of  the  downstream  half  of  the  plug. 
The  film  cooling  air  is  induced  through  three  sets  of  ports 
located  as  illustrated  in  Figure  5.  The  dimensional 
configuration  of  the  ports  is  shown  in  Figure  6.  The  second 
means  of  cooling  the  plug  provides  an  air  flow  to  convectively 
cool  the  downstream  half  of  the  plug.  This  is  accomplished 
with  a  shroud  that  directs  a  flow  of  air  along  the  backside 
of  the  plug  before  it  is  pulled  through  the  ports  to  be  used 
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for  film  cooling.  This  flow  of  air  is  called  tertiary  air 
flow  and  is  illustrated  in  Figure  5. 

B.  CUSP  MODIFICATION  OF  THE  PLUG 

The  symmetric  design  of  the  plug  was  modified  by  the 
addition  of  a  cusp  on  the  exhaust  gas  inlet  side  in  order  to 
provide  a  more  uniform  and  gradual  direction  of  the  gas  flow 
about  the  plug  and  to  physically  separate  the  flow  into  two 
symmetric  parts.  The  cusp  addition  is  picture  in  Figures  7 
and  9  and  dimensionally  illustrated  in  Figure  8.  The  cusp 
was  manufactured  from  0.102  cm  (0.04  in)  thick  sheet  aluminum. 
Each  side  of  the  cusp  was  formed  to  a  concave  shape  with  a 
24.77  cm  (9.75  in)  radius  of  curvature.  Axial  supports  were 
soldered  on  the  inside  of  the  cusp  to  prevent  warpage.  The 
cusp  was  then  glued  over  the  inlet  end  of  the  plug.  The 
apex  of  the  cusp  was  aligned  on  the  cross-sectional  center- 
line  of  the  mixing  stack  to  maintain  symmetry  of  the  flow 
about  each  half  of  the  mixing  stack.  The  apex  of  the  cusp 
extended  upstream  toward  the  entrance  of  mixing  stack  7.11 
cm  (2.8  in)  along  the  axial  centerline  from  the  end  of  the 
plug.  The  lip  formed  where  the  edges  of  the  cusp  met  the 
plug  was  filled  in  with  putty  to  provide  a  smooth  surface. 

The  corners  between  the  cusp  and  the  sides  of  the  mixing 
stack  were  sealed  with  silicone  epoxy  to  minimize  leakage. 


C.  MIXING  STACK  SIDE  MODIFICATION 


The  second  geometric  modification  to  Shaw’s  design  changed 
the  shape  of  the  curved  sides  of  the  mixing  stack.  In  Shaw's 
original  design,  at  the  position  where  the  two  air  flows 
around  the  plug  reconverged,  the  sides  of  the  mixing  stack 
had  a  137  degree  angle  bend  into  the  final  exit  duct  portion 
(see  Figure  5).  This  point  was  48.59  cm  (19.13  in)  from  the 
entrance  of  the  stack.  This  sharp  corner  was  removed  from 
the  stack  by  installing  a  curved  side  beginning  at  a  point 
43.18  cm  (17.0  in)  from  the  stack  entrance.  The  radius  of 
curvature  of  the  new  piece  was  43.18  cm  (17  in).  The  new 
sides  were  also  made  of  .102  cm  (.04  in)  sheet  aluminum. 

Edges  formed  where  the  new  section  was  fitted  into  the 
mixing  stack  were  filled  in  with  putty  to  maintain  a  smooth 
surface.  The  mixing  stack  with  the  side  modification  is 
dimensionally  shown  in  Figure  10.  The  picture  in  Figure  11 
highlights  the  area  changed  by  the  side  modification  and 
cusp  addition. 
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V.  EXPERIMENTAL  FACILITY 


Air  is  supplied  to  the  primary  nozzles  by  means  of  a 
centrifugal  compressor  and  associated  ducting  schematically 
illustrated  in  Figure  1.  The  mixing  stack  configuration 
being  tested  is  placed  inside  an  air  plenum  containing  an 
airtight  partition  so  that  two  separate  air  flows,  secondary 
and  tertiary,  may  be  measured.  The  air  plenum  facilitates 
the  accurate  measurement  of  secondary  and  tertiary  air  flows 
by  using  ASME  long  radius  flow  nozzles. 

A.  PRIMARY  AIR  SYSTEM 

The  circled  numbers  found  in  this  section  refer  to  lo¬ 
cations  on  Figure  1.  The  primary  air  ducting  is  constructed 
of  16-gage  steel  with  0.635  cm  (0.25  in)  thick  steel  flanges. 
The  ducting  sections  were  assembled  using  0.635  cm  (0.25  in) 
bolts  with  air  drying  silicone  rubber  seals  between  the 
flanges  of  adjacent  sections.  Entrance  to  the  inlet  ducting 
(T)  is  from  the  exterior  of  the  building  through  a  91.44  cm 
(3.0  ft)  square  to  a  30.48  cm  (1.0  ft)  square  reducer,  each 
side  of  which  has  the  curvature  of  a  quarter  ellipse.  A 
transition  section  T"  then  changes  the  30.48  cm  (1.0  ft) 
square  section  to  a  35.31  cm  (13.90  in)  diameter  circular 
section  ^  .  This  circular  section  runs  approximately 

9.14  m  (30  ft)  to  the  centrefugal  compressor  inlet. 


A  standard  ASME  square  edged  orifice  (4,  is  located  15 
diameters  downstream  of  the  entrance  reducer  and  11  diameters 
upstream  of  the  centrifugal  compressor  inlet,  thus  insuring 
stability  of  flow  at  both  the  orifice  and  compressor  inlet. 
Piezometer  rings  (?)  are  located  one  diameter  upstream  and 
one-half  diameter  downstream  of  the  orifice.  The  duct  sec¬ 
tion  also  contains  a  thermocouple  just  downstream  of  the 
orifice.  Primary  flow  is  measured  by  means  of  the  standard 
ASME  square  edged  orifice  designed  to  the  specifications 
given  in  the  ASME  power  test  code  [Ref.  5].  The  17.55  cm 
(6.902  in)  diameter  orifice  used  was  constructed  out  of  304 
stainless  steel  0.635  cm  (0.25  in)  thick.  The  inside  dia¬ 
meter  of  the  duct  at  the  orifice  is  35.31  cm  (13.90  in)  which 
yields  a  beta  (S  *  d/D)  of  0.497.  The  orifice  diameter  was 
chosen  to  give  the  best  performance  in  regard  to  pressure 
drop  and  pressure  loss  across  the  orifice  for  the  primary 
air  flow  rate  used  (1.71  Kg/sec  (3.77  lbm/sec)). 

The  centrifugal  compressor  7  used  to  provide  primary 
air  to  the  system  is  a  Spencer  Turbo  Compressor,  catalogue 
number  2S100-H,  rated  at  6000  cfm  at  2.5  psi  back  pressure. 
The  compressor  is  driven  by  a  three  phase,  440  volt,  100 
horsepower  motor. 

A  manually  operated  sliding  plate  variable  orifice  (6 
was  designed  to  constrict  the  flow  symmetrically  and  facili¬ 
tate  fine  control  of  the  primary  air  flow.  During  operation, 
the  butterfly  valve  (8s  ,  located  at  the  compressor's 
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discharge,  provided  adequate  regulation  of  primary  air  flow, 
eliminating  the  necessity  of  using  the  sliding  plate  valve. 
The  sliding  plate  valve  was  positioned  in  the  wide-open  po¬ 
sition  for  all  data  runs. 

On  the  compressor  discharge  side,  immediately  downstream 
of  the  butterfly  valve,  is  a  round  to  square  transition 
followed  by  a  90  degree  elbow  ^fo)  and  a  straight  section  of 
duct.  All  ducting  to  this  point  is  considered  part  of  the 
fixed  primary  air  supply  system.  A  transition  section 
is  fitted  to  this  last  square  section  which  reduces  the  duct 
cross  section  to  a  circular  section  29.72  cm  (11.17  in)  in 
diameter.  This  circular  ducting  tapers  down  to  a  diameter 
of  26.30  cm  (11.5  in)  to  provide  the  primary  air  inlet  to 
the  eductor  system  being  tested.  The  transition  is  located 
far  enough  upstream  of  the  model  to  insure  that  the  flow- 
reaching  the  model  is  fully  developed. 

B.  SECONDARY  AIR  PLENUM 

The  secondary  air  plenum,  pictured  in  Figure  1  and  2,  is 
constructed  of  1.905  cm  (0.75  in)  plywood  and  measures  1.22 
m  by  1.22  m  by  1.88  m  (4  ft  by  4  ft  by  6.17  ft).  It  serves 
as  an  enclosure  that  can  contain  all  or  only  part  of  the 
eductor  model  and  still  allow  the  exit  plane  of  the  mixing 
stack  to  protrude.  The  purpose  of  the  secondary  air  plenum 
is  to  serve  as  a  boundary  through  which  secondary  air  for 
the  eductor  system  must  flow.  Long  radius  ASME  flow  nozzles, 
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designed  in  accordance  with  ASME  power  test  codes  [Ref.  5] 
and  constructed  of  fiberglass,  penetrate  the  secondary  air 
plenum,  thereby  providing  the  sole  means  for  metering  the 
secondary  ail  reaching  the  eductor.  Appendix  D  of  Reference 
[1]  outlines  the  design  and  construction  of  the  secondary 
air  flow  nozzles.  By  measuring  the  temperature  of  the  air 
entering  and  the  pressure  differential  across  the  ASME  flow 
nozzles,  the  mass  flow  rate  of  secondary  air  can  be  deter¬ 
mined.  Flexibility  is  provided  in  measurement  of  the  mass 
flow  rate  of  secondary  air  by  employing  flow  nozzles  with 
three  different  throat  diameters:  20.32  cm  (8  in),  10.16  cm 
(4  in),  and  5.08  cm  (2  in).  By  using  a  combination  of  flow 
nozzles,  a  wide  variety  of  secondary  cross  sectional  areas 
can  be  obtained. 

A  secondary  air  flow  straightener ,  shown  in  Figures  1  and 
2,  consisting  of  a  double  screen  is  installed  1.22  m  (4  ft) 
from  the  open  end  of  the  secondary  air  plenum,  between  the 
ASME  long  radius  nozzles  and  the  primary  air  flow  nozzles. 

The  purpose  of  the  straightener  is  to  reduce  any  swirl  effect 
that  could  result  when  only  a  small  secondary  air  flow  area 
exists. 

C.  TERTIARY  AIR  PLENUM 

The  tertiary  air  plenum,  pictured  in  Figures  12  and  13, 
is  constructed  of  1.90  cm  (0.75  in)  plywood  and  measure  1.22 
m  by  1.22  m  by  1.22  m  (4ft  by  4  ft  by  4  ft).  It  serves  as 
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an  enclosure  that  completely  surrounds  the  mixing  stack  and 
allows  the  exit  and  entrance  regions  to  protrude.  An  air¬ 
tight  rubber  diaphragm  type  seal,  schematically  illustrated 
in  Figure  2  and  pictured  in  Figures  14  and  15,  is  located 
at  each  end  of  the  enclosure.  This  allows  measurement  of  a 
tertiary  air  flow  independent  of  the  secondary  air  flow. 
Tertiary  air  flow  is  measured  with  the  use  of  long  radius 
ASME  flow  nozzles  designed  in  accordance  with  ASME  test  codes 
[Ref.  5]  and  constructed  of  fiberglass.  These  nozzles  are 
located  so  that  they  penetrate  the  airtight  tertiary  air 
plenum,  thereby  providing  the  sole  means  for  metering  the 
tertiary  air  reaching  the  eductor.  By  measuring  the  tempera¬ 
ture  of  the  air  entering  and  the  pressure  differential  across 
the  ASME  flow  nozzles,  the  mass  flow  rate  of  tertiary  air 
can  easily  be  obtained.  Flexibility  in  measuring  the  ter¬ 
tiary  flow  is  provided  by  employing  different  size  flow 
nozzles:  two  of  20.32  cm  (8  in)  throat  diameter,  three  of 

10.16  cm  (4  in)  throat  diameter,  and  two  of  5.08  cm  (2  in) 
throat  diamter.  By  using  various  combinations  of  these  flow 
nozzles,  a  wide  variety  of  tertiary  cross  section  flow  areas 
can  be  obtained. 

The  interior  of  the  tertiary  air  plenum  is  pictured  in 
Figure  13.  The  stand  which  holds  the  mixing  stack  can  be 
seen  mounted  inside  the  plenum.  Figure  15  shows  the  diaphragm 
air  seal  at  the  exit  plane  of  the  mixing  stack.  As  can  be 
seen,  removalble  sections  were  located  in  the  exit  plane 
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door  to  allow  for  adjustments  to  the  mixing  stack  and  instru¬ 
mentation  without  removing  the  diaphragms. 

D.  ALIGNMENT 

The  alignment  of  the  mixing  stack  with  the  primary  air 
flow  nozzles  was  accomplished  using  a  level,  a  30.78  (12.0 

in)  rule  graduated  in  0.25  mm  (0.01  in)  and  45.72  x  30.48  cm 
(18.0  x  12.0  in)  square.  The  graduated  rule  and  square  were 
to  establish  the  stand-off  distance  (D/S)  and  to  center  the 
primary  flow  nozzles  within  the  entrance  area  of  the  mixing 
stack.  The  geometric  alignment  was  checked  for  accuracy 
using  pressure  readings  at  symmetric  points  on  the  model. 
Additional  verification  was  obtained  by  subsequent  symmetric 
exit  velocity  profiles  measurements.  The  three  axis  mounting 
stand,  pictured  in  Figure  16  allowed  alignment  adjustments 
to  be  performed  easily. 

E.  INSTRUMENTATION 

Pressure  taps  for  measuring  gage  pressures  are  located 
inside  the  primary  air  uptakes  just  prior  to  the  primary 
nozzles,  inside  the  secondary  air  plenum,  inside  the  ter¬ 
tiary  air  plenum,  and  at  various  points  on  the  model.  A 
variety  of  manometers,  pictured  in  Figure  17,  were  used  to 
indicate  the  pressure  differentials.  A  schematic  represen¬ 
tation  of  the  pressure  measuring  instrumentaion  is  illustrated 
in  Figures  18  and  19.  Monitoring  of  each  of  the  various 
pressures  was  facilitated  by  the  use  of  a  scanivalve  and  a 
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multiple  valve  manifold.  The  scanivalve  was  used  to  select 
the  pressure  tap  to  be  read,  while  the  multiple  valve  mani¬ 
fold  allowed  selection  of  the  optimum  manometer  for  the 
pressure  being  recorded.  A  vent  was  included  in  the  multiple 
valve  manifold  which  provided  a  means  of  venting  the  mano¬ 
meters  between  pressure  readings.  When  taking  readings  of 
the  pressure  distribution  in  the  mixing  stack,  it  was 
necessary  to  manually  change  the  tubing  from  one  end  of  the 
manometer  to  the  other  in  order  to  get  the  negative  pressure 
readings.  The  valve  manifold  provided  a  selection  of  a 
15.24  cm  (6.0  in)  inclined  water  manometer,  a  5.08  cm  (2.0 
in)  inclined  water  manometer,  and  a  1.27  cm  (0.5  in)  inclined 
oil  manometer  (specific  gravity  0.827).  In  addition,  the 
following  dedicated  manometers  were  used  in  the  system:  a 

50.80  cm  (20  in)  single  column  water  manometer  connected  to 
the  primary  air  flow  just  prior  to  the  primary  nozzles,  a 
1.27  cm  (50  in)  U-tube  water  manometer  with  each  leg  connec¬ 
ted  to  a  piezometric  ring  on  either  side  of  the  orifice 
plate  in  the  air  inlet  duct,  and  a  2.54  cm  (1.0  in)  inclined 
water  manometer  connected  to  the  upstream  piezometric  ring. 

Primary  air  temperatures,  measured  at  the  orifice  outlet 
and  just  prior  to  the  primary  nozzles,  are  measured  with 
copper-constantan  thermocouples.  The  thermocouples  are  in 
assemblies  manufactured  by  Honeywell  under  the  trade  name 
Megapak.  Polyvinyl  covered  20  gage  copper-constantan 
extension  wire  is  used  to  connect  the  thermocouples  to  a 


Newport  Digital  Pyrometer,  model  number  267,  which  provides 
a  digital  display  of  the  measured  temperature  in  degrees 
Fahrenheit.  The  Newport  Digital  Pyrometer  failed  in  use  and 
was  replaced  by  an  Omega  Digital  Thermometer,  Model  Number 
2176A.  Secondary/tertiary  ambient  air  temperature  is  mea¬ 
sured  with  a  mercury-glass  thermometer  and  recorded  in 
degrees  Fahrenheit. 

Velocity  profiles  at  the  mixing  stack  exit  plane  are 
obtained  by  using  a  pitot  tube,  pictured  in  Figure  20.  The 
tube  is  affixed  to  a  mounting  template  which  allows  accurate 
determination  of  the  major  axis,  minor  axis  ,  minor  and 
diagonal  positions  and  distances.  Alignment  pins  allow 
changes  in  velocity  traverse  directions  to  be  made  fast  and 
accurately.  The  pitot  tube  is  used  in  conjunction  with  the 
15.24  cm  (6  in)  inclined  water  manometer  for  obtaining  the 
velocity  pressure  head. 

An  additional  objective  of  this  research  was  to  assemble 
an  automated  data  acquisition  system.  The  data  which  was 
used  in  this  thesis  was  taken  "manually",  i.e.,  pressure 
and  temperature  values  visually  acquired  from  manometer  and 
thermometer.  Each  piece  of  data  was  then  inputed  into  a 
computer  for  reduction.  Although  equipment  limitation  pre¬ 
vented  full  automation,  "semi-automatic"  data  acquisition 
now  exists  and  is  described  in  Appendix  B. 


VI.  EXPERIMENTAL  METHOD 


Evaluation  of  the  eductor  model  requires  the  experimental 
determination  of  pressure  differentials  across  the  ASME  long 
radius  flow  nozzles,  temperature  of  primary  and  induced  air 
flows,  internal  mixing  stack  pressures,  and  mixing  stack  exit 
velocities.  These  experimentally  determined  quantities  are 
then  reduced  to  obtain  pumping  coefficients,  induced  air 
flow  rates,  pressure  distributions  within  the  mixing  stack, 
and  mixing  stack  exit  velocity  profiles.  In  addition  an 
indication  of  flow  behavior  in  the  stack  was  accomplished  by 
flow  visualization.  The  performance  characteristics  of  the 
eductor  model  are  then  evaluated  to  determine  the  model's 
relative  effectiveness. 

The  following  discussion  addresses  the  individual  perfor¬ 
mance  characteristics  of  the  eductor  model  and  how  they  were 
determined. 

A.  PUMPING  COEFFICIENTS 

The  secondary  pumping  coefficient  and  the  tertiary  pump¬ 
ing  coefficient  provide  the  basis  for  analysis  of  the  eductor 
model's  pumping  performance.  Thus,  changes  in  eductor  model 
parameters  which  affect  pumping  can  be  noted  by  a  change  in 
pumping  coefficient.  The  pumping  coef f icient (s)  is  desired 
at  the  operating  point  corresponding  to  that  of  the  ship¬ 
board  prototype.  At  this  point  there  is  no  restriction  of 
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the  secondary  (or  tertiary)  air  flow.  In  the  model,  this  is 
simulated  by  completely  opening  the  air  plenum(s)  to  the 
environment.  However,  at  this  condition,  the  secondary  (or 
tertiary)  air  flow  rate  cannot  be  measured.  Therefore,  the 
eductor's  characteristics  are  determined,  plotted,  and  then 
extrapolated  to  find  the  operating  point,  see  Figure  21. 

The  pumping  characteristics  of  the  eductor  model  are 
established  by  varying  the  associated  induced  air  flow  rate, 
either  secondary  or  tertiary,  from  zero  to  its  maximum 
measurable  rate.  This  rate  is  determined  by  sequentially 
opening  the  ASME  flow  nozzles  mounted  in  the  appropriate 
plenum  and  recording  the  pressure  drop  accross  the  nozzles. 
Values  for  nozzle  cross  sectional  area,  pressure  drop,  and 
induced  air  temperature  are  then  used  to  calculate  the 

0  .  44 

dimensionless  parameters  p*/T*  and  W*T*  or  PT*/TT*  and 

0-44 

WT*TT*  as  described  in  Appendix  A.  The  dimensionless 
parameters  are  then  plotted  as  illustrated  in  Figure  21. 
Extrapolation  of  the  pumping  characteristics  curve  to  inter¬ 
sect  with  the  zero  pressure/temperature  coefficient  abscissa 
locates  the  appropriate  operating  point  coefficient  of  the 
model . 

B.  INDUCED  AIR  FLOWS 

Two  induced  air  flows  are  identified  in  this  study:  se¬ 
condary  and  tertiary. 
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The  secondary  air  flow  is  the  amount  of  air  induced  by 
the  primary  nozzles  which  mixed  with  the  primary  air  flow 
to  reduce  the  exhaust  gas  temperature. 

The  tertiary  air  flow  indicates  the  amount  of  air  induced 
by  the  stack  pressure  distribution  to  provide  cooling  air 
through  the  ports  in  the  stack  plug.  In  addition  to  the 
method  of  metering  tertiary  air  flow  as  described  in  Section 
III,  the  effect  on  the  tertiary  pumping  coefficient  by  vary¬ 
ing  the  number  of  rows  of  cooling  ports  open,  see  Figure  6, 
was  investigated. 

C.  PRESSURE  DISTRIBUTIONS 

The  mixing  stack  axial  static  pressure  was  obtained  using 
a  series  of  pressure  taps  fixed  to  the  mixing  stack.  These 
taps  were  placed  in  two  axial  rows,  the  rows  being  along  the 
centerline  of  the  flow.  Along  each  row  the  taps  were  axially 
spaced  in  increments  of  one  quarter  primary  dimension.  The 
exact  location  of  the  pressure  taps  is  indicated  on  Figure 
2  and  4.  The  non-dimensional  stack  pressure,  PMS*,  (see 
Appendix  A)  is  plotted  versus  X/S  to  obtain  a  mixing  stack 
pressure  distribution,  see  Figures  26a  through  261. 

D.  EXIT  VELOCITY  PROFILES 

Velocity  profiles  at  the  mixing  stack  exit  were  calcula¬ 


ted  from  the  pressures  measured  using  the  pitot  tube  pictured 
in  Figure  20.  Since  it  was  impractical  to  obtain  a  complete 
three-dimensional  plot  of  velocities  at  the  exit  plane  of 


the  mixing  stack,  advantage  was  taken  of  the  symmetry  of 
the  velocity  surface  resulting  from  the  arrangement  of  the 
primary  nozzles.  Only  three  traverses  were  made.  The  first 
traverse  passes  vertically  along  the  centerline  down  the 
mixing  stack  exit,  the  second  traverse  passes  horizontally 
across  the  mixing  stack.  This  traverse,  along  the  center- 
line,  passes  across  the  top  of  the  plug.  The  third  traverse 
passes  diagonally  through  the  intersection  of  the  centerlines. 
Figure  22  illustrates  the  orientation  and  identification  of 
the  three  velocity  traverses,  while  Figures  27a  through  i, 
show  plots  of  the  velocity  traverses. 

E.  FLOW  VISUALIZATION 

Based  upon  Shaw's  results  it  was  felt  necessary  to 
obtain  a  visual  indication  of  the  flow  patterns  in  the  mixing 
stack.  Visualization  of  flow  was  accomplished  by  using 
string  or  yarn  placed  in  the  flow  stream.  Two  procedures 
were  used.  The  first  utilized  a  string  the  length  of  the 
mixing  stack.  With  this  string  the  shape  of  the  air  stream 
flow  throughout  the  entire  mixing  stack  could  be  seen.  The 
second  procedure  involved  the  insertion  of  a  probe  with  a 
small  string  attached  to  the  end  of  it  into  the  mixing  stack. 
With  this  procedure  it  was  very  easy  to  investigate  specific 
location  within  the  stack  for  flow  action,  particularly  for 
transition  regions  from  positive  to  negative  pressure.  In 
order  to  record  the  flow  visualizations  it  was  necessary  to 
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use  a  16  nun  motion  picture  camera  set  for  a  speed  of  60 
frames  per  second. 
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VII.  DISCUSSION  OF  EXPERIMENTAL  RESULTS 


Exhaust  eductor  systems  designed  for  marine  gas  turbine 
applications  must  substantially  cool  exhaust  gases,  present 
an  exterior  stack  surface  temperature  which  will  not  give 
an  easily  detectable  infrared  signature,  and  effectively  dis¬ 
burse  exhaust  gases.  In  order  to  quantitively  evaluate  the 
overall  eductor  model  performance,  four  areas  of  performance 
were  identified:  the  amount  of  secondary  air  flow  induced 
by  the  primary  air  flow,  referred  to  here  as  pumping;  the 
degree  of  mixing  of  primary  and  induced  air  flow  within  the 
mixing  stack  system,  referred  to  here  as  mixing;  the  amount 
of  uptake  back  pressure  impressed  upon  the  turbine  exhaust 
by  the  eductor  system;  and  the  amount  of  cooling  air  available 
to  reduce  the  exterior  stack  temperature  of  the  eductor 
system.  A  qualitative  means  of  assessing  flow  patterns  in 
the  mixing  stack  was  accomplished  by  flow  visualization 
with  string  placed  in  the  gas  flow  regions. 

The  eductor  model  and  its  modified  configurations  in 
this  study  were  designed  to  shield  the  primary  flow  nozzles 
from  an  overhead  view  and  to  provide  cooling  to  those  sur¬ 
faces  which  are  visible  from  above. 

The  initial  tests  performed  were  for  flow  visualization. 
The  tests  showed  the  following  condition  within  the  mixing 
stack.  First  there  appeared  to  be  large  scale  turbulence 
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within  the  stack.  This  was  demonstrated  by  the  wide  range 
of  travel  of  a  string  running  the  length  of  the  stack. 
Secondly,  there  was  evidence  of  some  sort  of  reverse  flow, 
probably  a  vortex  formation  on  the  downstream  sides  of  the 
plug.  Finally,  there  was  a  demonstration  of  flow  oscillatory 
interaction  at  the  outlet  end  of  the  plug.  The  range  of 
motions  of  the  strings  in  the  flow  is  shown  in  Figure  23. 

Further  model  testing  was  conducted  using  an  eductor 
model  with  the  addition  of  a  cusp  on  the  inlet  end  of  the 
symmetric  plug.  The  tests  were  conducted  with  a  primary 
flow  nozzle  area  to  mixing  stack  area  ratio  of  3.0,  a  mixing 
stack  length  L/S  of  3.0  and  a  standoff  distance  D/S  of  0.5. 
These  ratios  of  mixing  stack  length,  standoff  distance  and 
primary  nozzle  to  mixing  stack  area  were  used  as  a  result 
of  the  previous  study  by  Shaw  [Ref.  8].  The  reults  of  this 
first  set  of  eductor  system  tests  are  shown  in  Tables  I 
through  XIII. 

These  tests  showed  the  following  results.  First,  the 
uptake  back  pressure,  over  the  variation  in  secondary  and 
tertiary  flows  tested,  showed  a  reduction  of  up  to  .5  inches 
of  H20.  However,  at  the  operating  point  of  the  eductor 
system  (secondary  and  tertiary  areas  wide  open)  the  back 
pressure  was  the  same  as  Shaw's  9.85  inches  of  H20.  This 
uptake  back  pressure  may  be  too  high  and  would  impact  on 
the  turbine  efficiency.  Secondly,  the  eductor  model  pumping 
coefficient  was  improved  by  approximately  sixteen  percent 


over  that  of  Shaw.  Figure  24a  demonstrates  that  with  Shaw's 
design  the  pumping  coefficient  at  the  operating  point  would 
be  less  than  0.3.  With  the  addition  of  the  cusp  alone,  the 
pumping  coefficient  is  increased  above  0.35.  Finally,  the 
velocity  profiles  in  all  three  direction;  vertical,  horizon¬ 
tal,  and  diagonal  showed  little  change  from  Shaw's  data. 

On  the  average  these  profiles  were  relatively  flat.  These 
velocity  traverses  are  depicted  in  Figures  17a  through  2 7 i . 

The  next  mixing  stack  modification  was  the  replacement 
of  a  sharp  corner  in  the  sides  with  a  curved  section  that 
smoothly  transitioned  from  the  ducting  over  the  plug  to  the 
exit  ducting.  The  results  of  the  tests  performed  after  this 
modification  are  tabulated  in  Tables  XIII  through  XVIII. 

Initial  comparisons  of  data  were  with  the  back  pressures 
obtained.  Over  the  ranges  of  secondary  and  tertiary  flow 
combinations  the  back  pressures  increased  over  those  values 
recorded  for  the  cusp  only  addition.  However,  at  the  wide 
open  operating  point  the  back  pressure  was  the  same  as  that 
with  the  cusp  only  and  with  Shaw's  data,  9.85  inches  of  H20. 
With  the  addition  of  the  side  modification  the  pumping  co¬ 
efficient,  over  the  range  of  flows  used,  was  reduced  from 
those  values  obtained  without  the  side  modification  (see 
Figure  24c) .  The  pumping  coefficients  at  the  operating 
point  without  and  with  the  side  modification  were  essentially 
the  same.  Finally,  the  velocity  profiles  in  all  three 
directions  showed  more  variation  with  the  side  modification 
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than  previously.  However,  in  the  aggragate  these  profiles 
were  still  relatively  flat.  These  traverses  are  in  Figures 
27g,  h  and  i. 

The  following  observations  were  made  in  comparing  the 

0*44 

cooling  (tertiary)  pumping  coefficients  (WT*TT*  ).  Con¬ 
comitant  with  an  increase  or  decrease  in  secondary  pumping 
coefficient  was  an  inverse  action  with  the  tertiary  pumping 
coefficient.  The  addition  of  the  cusp  decreased  the  tertiary 
pumping  coefficient  from  .060,  the  value  obtained  by  Shaw, 
by  approximately  33%  to  .040  (see  Figure  25c).  Over  the 
ranges  of  flow  combinations  tested  the  addition  of  the  side 
modification  improved  the  tertiary  pumping  coefficient.  At 
the  operating  point  the  values  of  the  tertiary  pumping  co¬ 
efficients  for  with  and  without  the  side  modification 
essentailly  converged.  These  results  are  displayed  in  Fi- 
fures  25a  through  25g. 

It  was  further  noted  that  decreasing  the  number  of 
tertiary  ports  open  had  the  effect  of  increasing  the 
secondary  pumping  coefficient  and  consequently  decreasing 
the  tertiary  pumping  coefficient. 

The  final  indicator  looked  at  was  the  pressure 
distribution  in  the  mixing  stack.  The  data  obtained  was 
essentially  unchanged  from  that  of  Shaw's.  Positive  pressures 
occurred  along  the  exterior  curved  surface  of  the  mixing 
stack.  As  with  Shaw,  a  large  negative  pressure  (PMS*)  = 

-  0.30,  occurs  at  the  junction  of  the  curved  surface  and 


the  rectangular  mixing  stack  exit  plane.  It  was  noted  that 

with  the  addition  of  the  cusp  only,  this  large  negative 

pressure  (PMS*)  increased  to  a  value  of  -0.35  (Figure  26a). 

With  the  further  modification  of  the  curved  side  the  large 

negative  pressure  was  significantly  reduced  to  -.15,  a  value 

( 

much  less  than  Shaw's  (Figure  26i).  No  significant  pressures 
were  noted  along  the  flat  surface  of  the  mixi’ng  stack. 


VI.  CONCLUSIONS 


This  investigation  studied  the  effects  on  the  eductor 
system’s  overall  performance  of  a  mixing  stack  geometry 
which  has  a  rectangular  cross-section  and  employs  a  plug 
placed  over  the  primary  flow  nozzles.  The  plug  was  modified 
with  the  addition  of  a  cusp  on  the  inlet  end.  The  sides  of 
the  mixing  stack  was  further  modified  to  remove  a  sharp 
corner. 

The  conclusions  resulting  from  this  investigation  are 
given  here. 

1.  The  addition  of  the  cusp  on  the  plug  increased  the 
secondary  pumping  coefficient  and  increased  a  ne¬ 
gative  pressure  at  a  point  within  the  mixing  stack. 

2.  The  addition  of  the  cusp  reduced  the  tertiary  cool¬ 
ing  pumping  coefficient. 

3.  The  modification  of  the  sides  of  the  mixing  stack 
reduced  the  negative  pressure  at  the  specific  point 
within  the  mixing  stack. 

4.  The  modification  to  the  curved  side  of  the  mixing 
stack  did  not  improve  the  performance  of  the  eductor 
system  and  should  be  removed. 
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VII.  RECOMMENDATIONS  FOR  FURTHER  STUDY 


This  study  showed  the  effects  on  a  specific  eductor 
system  with  several  geometric  configurations.  This  research 
adds  some  information  to  the  overall  data  base  upon  which 
further  research  can  be  conducted.  Recommendations  for 
further  investigation  of  this  eductor  system  are  presented 
here . 

1.  Further  explore  the  flow  patterns  within  the  stack 
using  tufting  or  some  other  flow  visualization 
process . 

2.  Test  the  eductor  model  with  the  following  modifica¬ 
tions:  remove  the  side  modification,  add  an  exterior 
shroud  to  the  mixing  stack  and  add  a  diffusor  con¬ 
figuration  to  the  exit  plane  of  the  mixing  stack. 
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FIGURE  2.  EDUCTOR  MODEL  WITH  SECONDARY  AND  TERTIARY  PLENUMS 


FIGURE  3.  SCHEMATIC  OF  MIXING  STACK  EDUCTOR 


FIGURE  8.  DIMENSIONAL  ILLUSTRATION  OF  MIXING  STACK  WITH  CUSP 


CURE 


FIGURE  11.  MIXING  STACK  WITH  CUSP  AND  SIDE  MODIFICATION 
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FIGURE 


I7-  mixing  stack 
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FIGURE  18.  SCHEMATIC  OF  INSTRUMENTATION 


FIGURE  19.  SCHEMATIC  OF  INSTRUMENTATION  FOR  PRIMARY  AIR  FLOW  MEASUREMENT 
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EXPERIMENT I AL  PUMPING 
COEFFICIENT 


FIGURE  21.  EXAMPLE  PLOT  OF  PUMPING  COEFFICIENT 
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FIGURE  24a 


EXPERIMENTAL  PUMPING 
COEFFICIENT  COMPARISON 


>*,“0  -P  -P 
L  C  £.  L 

O  O  &  <0 

-P  -P 

-P  CL 

t_  a  cm  m 

a  3  i  i 

-p  o  *  + 


(W*)  (T*)  \  44 

FIGURE  24b 
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FIGURE  24o 


EXPERIMENTAL  FILM  PUMPING 
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<WT*>  (TT*)  \  44 
FIGURE  25a 


EXPERIMENTAL  FILM  PUMPING 
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(WT*>  (TT*)  44 

FIGURE  25b 


EXPERIMENTAL  FILM  PUMPING 
COEFFICIENT  COMPARISON 
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FIGURE  25o 


EXPERIMENTAL  FILM  PUMPING 
COEFFICIENT  COMPARISON 
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(WT*)  CTT*) \  44 
FIGURE  25d 


EXPERIMENTAL  FILM  PUMPING 


FIGURE  25. 


EXPERIMENTAL  FILM  PUMPING 
COEFFICIENT  COMPARISON 
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FIGURE  25f 


EXPERIMENTAL  FILM  PUMPING 
COEFFICIENT  COMPARISON 
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FIGURE  25g 
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FIGURE  26a 
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AXIAL  DISTANCE  (X/S) 
FIGURE  26b 


AXIAL  PRESSURE 
DISTRIBUTION  COMPARISON 
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FIGURE  26o 


AXIAL  PRESSURE 


AXIAL  DISTANCE  (X/S) 

FIGURE  26d 


AXIAL  PRESSURE 
DISTRIBUTION  COMPARISON 


AXIAL  DISTANCE  (X/S) 

FIGURE  26e 


AXIAL  PRESSURE 
DISTRIBUTION  COMPARISON 
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AXIAL  DISTANCE  (X/S) 
FIGURE  26h 


AXIAL  PRESSURE 
DISTRIBUTION  COMPARISON 
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AXIAL  DISTANCE  (X/S) 
FIGURE  26 i 


AXIAL  PRESSURE 
DISTRIBUTION  COMPARISON 
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AXIAL  DISTANCE  (X/S) 
FIGURE  26 j 


AXIAL  PRESSURE 
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iXIAL  DISTANCE  (X/S) 
FIGURE  26k 


AXIAL  PRESSURE 


AXIAL  DISTANCE  (X/S) 
FIGURE  261 
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DIAGONAL  DISTANCE  (IN) 
FIGURE  27a 


VERTICAL  VELOCITY  TRAVERSE 


VERTICAL  DISTANCE  (IN) 
FIGURE  27o 


DIAGONAL  VELOCITY  TRAVERSE 


DIAGONAL  DISTANCE  (IN) 
FIGURE  27 d 


HORIZONTAL  VELOCITY  TRAVERSE 


HORIZONTAL  DISTANCE  (IN) 
FIGURE  27a 


VERTICAL  VELOCITY  TRAVERSE 


VERTICAL  DISTANCE  (IN) 
FIGURE  27f 


DIAGONAL  VELOCITY  TRAVERSE 
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FIGURE  27q 
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HORIZONTAL  DISTANCE  (IN) 
FIGURE  27h 


VERTICAL  VELOCITY  TRAVERSE 


VERTICAL  DISTANCE  (IN) 
FIGURE  27 i 
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APPENDIX:  A 


ONE -DIMENSIONAL  ANALYSIS  OF  A  SIMPLE  EDUCTOR 

This  investigation,  being  an  extension  of  the  work  of 
Shaw  [Ref.  8],  and  others  Ref.  1,  2  and  7]  uses  the  same 
one-dimensional  analysis  of  a  simple  eductor  system.  Simi¬ 
larity  between  the  basic  geometry  tested  by  previous  resear 
chers  was  maintained  in  order  to  correlate  data.  The  demen 
sionless  parameters  controlling  the  flow  phenomena  used 
previously  were  also  used  in  this  investigation  along  with 
the  basic  means  of  data  analysis  and  presentation.  Dynamic 
similarity  was  maintained  by  using  Mach  number  similarity 
to  establish  the  model's  primary  flow  rate. 

Although  the  analysis  presented  here  is  for  an  eductor 
model  with  only  primary  and  secondary  air  flows,  it  should 
be  kept  in  mind  that  many  of  the  results  presented  are  for 
systems  with  primary,  secondary,  and  tertiary  air  flows. 
Systems  with  tertiary  and  film  or  wall  cooling  air  flows 
have-  been  non-dimensionalized  with  the  same  base  parameters 
as  the  secondary  air  flow  and  have  been  calculated  using 
the  same  one-dimensional  analysis.  This  allows  for  easy 
comparison  of  the  results.  Parameters  pertaining  to  the 
secondary  systems  are  subscripted  with  an  "s",  those 
relating  to  the  tertiary  box  are  subscripted  with  a  "t". 
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I 


A.  MODELING  TECHNIQUE 

Dynamic  similarity  between  the  models  tested  and  the 
actual  prototype  was  maintained  by  using  the  same  primary 
air  flow  Mach  number.  For  the  primary  air  flow  Mach  number 
used  (0.064),  and  based  on  the  average  flow  properties  within 
the  mixing  stack  and  the  hydraulic  diameter  of  the  mixing 
stack,  the  air  flow  through  the  eductor  system  is  turbulent 
(Re>105).  As  a  consequence  of  this,  momentum  exchange  is 
predominant  over  shear  interaction,  and  the  kinetic  and 
internal  energy  terms  are  more  influential  on  the  flow  than 
are  viscous  forces.  It  can  also  be  shown  that  the  Mach 
number  represents  the  ratio  of  kinetic  energy  of  a  flow  to 
its  internal  energy  and  is,  therefore,  a  more  significant 
parameter  than  the  Reynolds  number  in  describing  the  primary 
flow  through  the  uptakes. 

B.  ONE -DIMENSIONAL  ANALYSIS  OF  A  SIMPLE  EDUCTOR 

The  theoretical  analysis  of  an  eductor  may  be  approached 
in  two  ways.  One  method  attempts  to  analyze  the  details  of 
the  mixing  process  of  the  primary  and  secondary  air  streams 
as  it  takes  place  inside  the  mixing  stack.  This  requires  an 
interpretation  of  the  mixing  phenomenon  which,  when  applied 
to  a  multiple  nozzle  system,  becomes  extremely  complex.  The 
other  method,  which  was  chosen  here,  analyzes  the  overall 
performance  of  the  eductor  system  and  is  not  concerned  with 
the  actual  mixing  process.  The  one-dimensional  analysis 
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based  on  a  single  primary  nozzle  exhausting  into  a  mixing 
stack,  as  shown  in  Figure  A-l.  To  avoid  repetition  with 
previous  reports,  only  the  main  parameters  and  assumption 
will  be  represented  here.  A  complete  derivation  of  analysis 
used  can  be  found  in  References  [1]  and  [4].  The  one-dimen¬ 
sional  flow  analysis  of  the  simple  eductor  system  described 
depends  on  the  simultaneous  solution  of  the  continuity, 
momentum  and  energy  equations  coupled  with  the  equation  of 
state,  all  compatible  with  specific  boundary  conditions. 

The  idealizations  made  for  simplifying  the  analysis  are 
as  follows: 


1.  The  flow  is  steady  state  and  incompressible. 

2.  Adiabatic  flow  exists  throughout  the  eductor  with 
isentropic  flow  of  the  secondary  stream  from  the 
plenum  (at  section  0)  to  the  throat  or  entiance  of 
the  mixing  stack  (at  section  1)  and  irreversible 
adiabatic  mixing  of  the  primary  and  secondary  streams 
occurs  in  the  mixing  stack  (between  sections  1  and  2) . 

3.  The  static  pressure  across  the  flow  at  the  entrance 
and  exit  planes  of  the  mixing-tube  (at  sections  1 
and  2)  is  uniform. 

4.  At  the  mixing-stack  entrance  (section  1)  the  primary 


flow  velocity  Up  and  temperature  Tp  are  uniform 
across  the  primary  stream,  and  the  secondary  flow 


velocity  Ug  and  temperature  Tg  are  uniform  across 
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the  secondary  stream,  but  Up  does  not  equal  Ug ,  and 

T  does  not  equal  T  . 

P  H  s 

5.  Incomplete  mixing  of  the  primary  and  secondary  streams 

in  the  mixing  stack  is  accounted  for  by  the  use  of  a 

non-dimensional  momentum  correction  factor  K  which 

m 

relates  the  actual  momentum  rate  to  the  pseudo-rate 

based  on  the  bulk-average  velocity  and  density  and 

by  the  use  of  a  non-dimensional  kinetic  energy 

correction  factor  K  which  relates  the  actual  kinetic 

e 

energy  rate  to  the  pseudo-rate  based  on  the  bulk- 
average  velocity  and  density. 

6.  Both  gas  flows  behave  as  perfect  gases. 

7.  Flow  potential  energy  position  changes  are  negligible. 

8.  Pressure  changes  PsQ  to  Pgl  and  P.^  to  P&  are  small 
relative  to  the  static  pressure  so  that  the  gas 
density  is  essentially  dependent  upon  temperature 
(and  atmospheric  pressure) . 

9.  Wall  friction  in  the  mixing  stack  is  accounted  for 
with  the  convectional  pipe  friction  factor  term 
based  on  the  bulk-average  flow  velocity  Um  and  the 
mixing  stack  wall  area  A . 

The  following  parameter,  defined  here  for  clarity,  will 
be  used  in  the  following  development. 


area  ratio  of  primary  flow  area  to  mixing 
stack  cross  sectional  area 


143 


J!L  area  ratio  of  wall  friction  area  to  mixing 

Am  stack  cross  sectional  area 

k  momentum  correction  factor  for  primary  mix- 

P  ing 

K  momentum  correction  factor  for  mixed  flow 

m 

f  wall  friction  factor 

Base  on  the  continuity  equation,  the  conservation  of 

mass  principle  for  steady  flow  yields 

W  =  W  +  W  +  W.  (1) 

m  p  s  t 


where 

W  =  p  U  A 
P  p  p  P 

W  =  P  U  A 
s  s  s  s 

wt  =  PtutAt  ('la> 

W  »  p  U  A 
m  m  m  m 

All  of  the  above  velocity  and  density  terms,  with  the 

exception  of  p  and  U  ,  are  defined  without  ambiguity  by 
mm 

the  virtue  of  idealizations  (3)  and  (4)  above.  Combining 
equations  (1)  and  (la)  above,  the  bulk  average  velocity  at 
the  exit  plane  of  the  mixing  stack  becomes 


U  =  — 


W  +  +  W 


where  Am  is  fixed  by  the  geometric  configuration  and 


where  T  is  calculated  as  the  bulk  average  temperature  from 
m 

the  energy  equation  (9)  below.  The  momentum  equation  stems 
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from  Newton's  second  and  third  laws  of  motion  and  is  the 
conventional  force  and  momentum-rate  balance  in  fluid  mechanics . 


w  u  w  u  w.u.  w  u 

K  (-2—2.)  +  (-L-1)  +  C— |— =■)  +  P  A  =  K  (JOj)  t  PA 

rv  rr  vrry  K  CT  '  ..  m  v  nr  '  - 


P  8, 


™  g. 


Note  the  introduction  of  idealizations  (3)  and  (5).  To 
account  for  a  possible  non-uniform  velocity  profiles  across 
the  primary  nozzle  exit,  the  momentum  correction  factor  Kp 
is  introduced  here.  It  is  defined  in  a  manner  similar  to 
that  of  Km  and  by  idealization  (4)  ,  supported  by  work  con¬ 
ducted  by  Moss,  it  is  set  equal  to  unity.  Kp  is  carried 
through  this  analysis  only  to  illustrate  its  effect  on  the 
final  result.  The  momentum  correction  factor  for  the 
mixing  stack  exit  is  defined  by  the  relation 


Km  "  W  U  J 

mm  ij 


U  2  P  dA 
m  2 


where  U  is  evaluated  as  the  bulk-average  velocity  from 
equation  (lb) .  The  wall  skin  friction  force  F^r  can  be 
related  to  the  flow  stream  velocity  by 


(5) 


using  idealization  (9) .  As  a  reasonably  good  approximation 
for  turbulent  flow,  the  friction  factor  may  be  calculated 
from  the  Reynolds  number 
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f  =  0.046(Rem) 


0  .  2 


(6) 


Applying  the  conservation  of  energy  principle  to  the  steady 
flow  system  in  the  mixing  stack  between  the  entrance  and 
exit  planes. 


U 


W  (h  +  J2-)  +  W  (h  +  O-)  +  W  (h  +  tX-) 

P  P  Z  sv  s  2g~}  t1-  t  2g7 


lO 
> 

7 

c 


u. 


u  2 

=  W  (h  +  K  JL_) 
m  elg^T 


(7) 


neglecting  potential  energy  of  position  changes  (idealization 
7).  Note  the  introduction  of  the  kinetic  energy  correction 


factor  Ke>  which  is  defined  by  the  relation 


r^7 

m  m 


A 

r  m 


U  3P  dA 
2  2 


(8) 


It  may  be  demonstrated  that  for  the  purpose  of  evaluat¬ 
ing  the  mixed  mean  flow  temperature  Tm,  the  kinetic  energy 


terms  may  be  neglected  to  yield 


hm  = 


w  w 

h  +  —  h  +  X  h. 
p  W_  s  W_  t 


m 


m 


(9) 


where  T  =  $(10  only,  with  the  idealization  (6). 
m  m 


The  energy  equation  for  the  isentropic  flow  of  the 
secondary  air  from  the  plenum  to  the  entrance  of  the  mixing 
stack  may  be  shown  to  reduce  to 


-  P  U  2 
o  s  _  s 

9  S  ^ 


(10) 
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similarly,  the  energy  equation  for  the  tertiary  air  flow 
reduces  to 


The  forgoing  equations  may  be  combined  to  yield  the 
vacuum  produced  by  the  eductor  action  in  either  the  secondary 
or  tertiary  air  plenums.  For  the  secondary  air  plenum,  the 
vacuum  produced  is 


-  P 


os 


pr(K 

’c  m 


W  2 

_e _ 

p  Vp 


w 


A  p 
SHS 


,  A  W  2 
r  -I  m-v  m 

1  2  A  J  A  o 

s  mm 


(K. 


m 


-  A 

1  Ji'n 

2  A  JJ 

m 


(ID 


where  it  is  understood  that  A^  and  pg  apply  to  the  secondary 
flow  at  this  same  section,  and  A  and  apply  to  the  mixed 
flow  at  the  exit  of  the  mixing  stack  system.  P  is  atmos- 
pheric  pressure,  and  is  equal  to  the  pressure  at  the  exit 
of  the  mixing  stack.  A  is  the  area  of  the  inside  wall  of 

W 

the  mixing  stack. 

For  the  tertiary  air  plenum,  the  vacuum  produced  is 
(W  +WJ  2 

ppp' 


Pa  “  Pot  g  A  ^p  (A  P  +A  1  +  AP.^'2  X"-1 
&cmrppss  tt  t 


V  f 

■  A-V'1.  *  1st” 

mm  m 


(11a) 


where  the  primary  flow  now  consists  of  both  the  primary 
and  secondary  air  flows. 
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C.  NON-DIMENSIONAL  FORM  OF  THE  SIMPLE  EDUCTOR  EQUATION 


‘s  =  T2  ' 


P 


* 

t 


As 


a  flow  density  ratio  of  the  tertiary 
or  film  cooling  flow  to  primary  flows. 
(Note  that  since  the  fluids  are  con¬ 
sidered  perfect  gases, 


an  area  ratio  of  secondary  flow  area 
to  primary  flow  area 


A? 


an  area  ratio  of  tertiary  flow  area 
to  primary  flow  area 


With  these  non-dimensional  groupings,  equations  (11)  and 


(11a)  can  be  rewritten  in  dimensionless  form.  Since  both 
equations  follow  the  same  format,  only  the  results  for  the 


secondary  air  plenum  will  be  presented  here. 


*  2  S^p  -  ^  '  w*'Kp 


m 


m 


m 


*  B)) 

r  pm 


where 


r  A 

8  "  Km  +  2  . 

m 


This  mav  be  rewritten  as 


(12) 


Cl  +  C2W*(T  +  1)  +  c3w*2t* 


(13) 
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where 


C'  =  25E(ICp  -  ^  8), 

m  *  m 
A  2 

C2  =  *  2  (  t-2.)  3,  and 
An 

«=•  -  4^  ■ 

m  pm 

As  can  be  ssen  from  equation  (13)  , 
AP*  =  F (W* ,  T*) . 


The  additional  dimensionless  quantitites  listed  below 
were  used  to  correlate  the  static  pressure  distribution  down 
the  length  of  the  mixing  stack. 


PMS* 


PMS 


s 

0"T 


a  pressure  coefficient  which  compares 
the  pumping  head  PMS  for  the  secondary 


flow  to  the  driving  head 


of  the 


primary  flow,  where  PMS  *  static 
pressure  along  the  mixing  stack  length 
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APPENDIX:  B 


FORMULA 


Presented  here  are  the  fomulas  used  to  obain  the  primary 
and  secondary  mass  flow  rates.  According  to  the  ASME  Power 
Test  Code  (5) ,  the  general  equation  for  mass  flow  rate  appear¬ 
ing  in  equation  (a) 

W(lbm/sec)  =  (0.12705)  K  A  Y  Fo  (p  AP) °  5  (a) 

a 

may  be  used  with  flow  nozzles  and  square  edge  orifices 
provided  the  flow  is  subsonic.  In  the  above  equation,  K 
(dimensionless)  represents  the  flow  coefficient  for  the 

.  -0.5 

metering  device  and  is  defined  as  K  =  C(1  -B  )  where  C 

is  the  coefficient  of  discharge  and  B  is  the  ratio  of  throat 

to  inlet  diameters;  A(in2)  is  the  total  cross  sectional  area 

of  the  metering  device;  Y  (dimensionless)  is  the  expansion 

factor  for  the  flow;  F  (dimensionless)  is  the  area  thermal 

21 

expansion  factor;  p  (lbm/ft3)  is  the  flow  mass  density;  and 
AP  (inches  H20)  is  the  differential  pressure  across  the 
metering  device.  Each  of  these  quantitites  are  evaluated, 
according  to  the  guidelines  set  forth  in  Reference  [5],  for 
the  specific  type  of  flow  measuring  device  used. 

Using  a  square  edge  orifice  for  measurement  of  the  pri¬ 
mary  mass  flow  rate,  the  quantities  in  equation  (a)  are 
defined  as  follows: 
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1.  The  flow  coefficient  K  is  0.62  based  on  a  8  of  0.S02 
and  a  constant  coefficient  of  discharge  over  the 
range  of  flows  considered  of  0.60. 

2.  The  orifice  area  is  37.4145  in2. 

3.  Corresponding  to  the  range  of  pressure  ratios  en¬ 
countered  across  the  orifice,  the  expansion  factor 
Y  is  0.98. 

4.  Since  the  temperature  of  the  metered  air  is  nearly 
ambient  temperature,  thermal  expansion  factor  is 
essentially  1.0. 

5.  The  primary  air  mass  density  PQr  is  calculated  using 
the  perfect  gas  relationship  with  pressure  and  tem¬ 
perature  evaluated  upstream  of  the  orifice. 

Substituting  these  values  into  equation  (a)  yields 

W  (lbm/sec)  =  (2.8882)  (PorAPor)0*5  (b) 

The  secondary  mass  flow  rate  is  measured  using  long 
radius  flow  nozzles  for  which  case  the  quantities  in  equation 
(a)  becomes: 

1.  For  a  flow  nozzle  installed  in  a  plenum,  3  is  approxi¬ 
mately  zero  in  which  case  the  flow  coefficient  is 
approximately  equal  to  the  coefficient  of  discharge. 
For  the  range  of  secondary  flows  encountered,  the 
flow  coefficient  becomes  0.98. 

2.  A  is  the  sum  of  the  throat  areas  of  the  flow  nozzles 


3.  Since  the  pressure  ratios  across  the  flow  nozzles 
are  very  close  to  unity,  the  expansion  coefficient  Y 
is  1.0. 

4.  Since  the  temperature  of  the  metered  air  is  nearly 
ambient  temperature,  the  thermal  expansion  factor  is 
essentially  1.0. 

5.  The  secondary  air  mass  density  P  is  evaluated  using 
the  perfect  gas  relationship  at  ambient  conditions. 

Substituting  these  values  into  equation  (a)  yields  the 
equation  for  the  secondary  mass  flow  rate  measured  using 
long  radius  flow  nozzles. 

o .  s 

Ws  (lbm/sec)  =  (0.13451)  A  (PgAP,.)  (c) 
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APPENDIX:  C 


DATA  ACQUISITION  AUTOMATION 

The  method  used  in  taking  data  with  this  and  previous 
theses  has  been  completely  manual.  Although  the  procedure 
has  been  adequate,  it  has  been  combersome.  An  ancillary 
objective  of  this  thesis  has  been  to  assemble  a  system  to 
automatically  collect,  reduce  ans  store  the  data  required. 

The  system  developed  utilizes  the  following  equipment: 

1)  Hewlett  Packard  Model  8S  Computer 

2)  Hewlett  Packard  Model  82901M  Flexible  Disc 
Drive 

3)  Hewlett  Packard  Model  722S-A  Plotter 

4)  Hewlett  Packard  Model  3497A  Data  Acquisition/ 
Control  Unit. 

These  equipments  were  connected  via  a  Hewlett  Packard  Model 
82937A  Interface  Bus.  Additionally,  the  HP-85  had  installed 
the  Mass  Storages,  Plotter/Printer,  and  Matrix  ROMS  (Read 
Only  Memories)  .  Figures  C-I  and  C-II  show  a  system  schematic 
and  the  actual  equipment  used. 

Pressure  reading  were  taken  with  poly-flo  tubing  from 
the  scanivalve  via  a  manifold  arrangement  to  Validyne  Model 
DP  103  pressure  transducers.  Two  transducers  are  utilized, 
one  with  a  pressure  range  of  0  to  5.5  inches  of  H20  (0-0.2 
psi)  ,  the  other  with  a  pressure  range  of  0  to  35  inches  of 
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H20  (0-1.2  psi).  The  DC  electrical  signal  from  the  trans¬ 
ducer  is  amplified  through  a  Validyne  CD10  Carrier  Demodula¬ 
tor.  The  carrier  demodulator  also  provides  controls  for 
setting  the  zero  point  of  the  transducer  and  varying  the 
span  of  the  voltage  output  from  the  unit.  The  DC  output 
from  the  demodulator  is  measured  by  the  voltmeter  within  the 
HP-3497A  data  acquisition  unit.  The  3497A  transmits  the 
voltage  value  to  the  HP-85  where  a  subroutine  of  the  acqui¬ 
sition  program  converts  the  voltage  to  a  pressure  value. 

Temperature  readings  are  taken  in  a  similar  manner.  The 
thermocouple  voltages  are  read  directly  bu  the  349 7A  data 
acquisition  unit.  A  different  subroutine  within  the  data 
acquisition  program  converts  the  voltage  reading  to 
temperature  values. 

The  conversion  formula  for  temperature  was  derived  from 
reference  tables  for  type  T  copper-constantan  thermocouples. 
Voltages  from  the  reference  table  at  ten  degree  increments 
from  0  to  120  degrees  Fahrenheit  were  input  into  a  program 
from  the  HP-85  Math  Pac.  This  program  provided  the 
coefficients  for  a  degree  Chebyshev  polynomial.  The  following 
formula  results: 

T  =  .0608773(V) 3  -  1 . 289098 (V) 2  +  46.58608(V)  +  32.00443 
where  T  is  degrees  Fahrenheit  and  V  is  millivolts. 

Unlike  the  thermocouple,  calibration  of  the  transducers 
must  be  accomplished  everytime  the  equipment  is  energized. 

This  is  because  the  location  of  the  zero  point  is  not  always 
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constant.  Calibration  of  the  transducers  required  use  of 
the  building  air  supply.  The  compressor  in  the  building  is 
a  Worthington  Air  Compressor  set  at  a  nominal  operating 
pressure  of  100  psig.  Sufficient  air  pressure  was  obtained 
by  regulating  the  supply  valve  from  the  building  air  system 
and  installing  a  throttling  valve  prior  to  a  manifold.  Both 
the  pressure  transducers  and  a  20  inch  water  manometer  are 
connected  to  the  manifold.  In  addition  an  exhaust  valve  that 
can  be  throttled  is  also  attached  to  the  manifold.  With  the 
building  air  supply  completely  disconnected  and  all  manifold 
valves  open,  the  zero  point  is  set.  The  air  supply  is  then 
connected  to  the  manifold  and  the  pressure  is  raised  to  5.0 
or  20  inches  of  water  depending  on  the  transducer  to  be 
calibrated.  A  vacuum  can  be  obtained  by  inserting  a  tee 
fitting  into  the  poly-flo  tubing  ahead  of  the  inlet  throttling 
valve.  This  forms  an  eductor  that  takes  a  suction  on  the 
manifold.  By  varying  the  building  air  supply  valve  and  the 
throttle  valves  connected  to  the  manifold  the  desired  vacuum 
can  be  reached.  When  calibrating  the  0  to  35  inch  transducer 
it  is  vital  that  the  valve  isolating  the  lower  range  trans¬ 
ducer  be  secured.  An  additional  safety  precaution  is  to 
have  a  vent  valve  installed  between  the  manifold  and  the 
transducer.  This  safety  valve  is  to  be  open  whenever  the 
5.5  inch  transducer  is  not  in  use.  After  raising  the  pressure 
to  the  high  end  calibration  point  the  pressure  is  trapped  in 
the  manifold  by  closing  all  valves.  A  transducer  reading  is 
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then  taken  and  recorded  with  the  pressure  reading  from  the 
manometer.  The  process  is  then  repeated  with  several  values 
ath  the  low  end  of  the  pressure  range.  These  voltages  and 
corresponding  pressure  values  are  then  inputed  into  the  Che- 
byshev  polynomial  program  and  the  coefficients  for  a  1st 
degree  euqation  are  obtained.  These  coefficients  are  then 
inserted  into  the  data  acquisition  program. 

It  is  important  to  note  that  the  transducer  can  be  set 
for  reading  either  positive  or  negative  differential  pressure. 
Once  a  transducer  has  been  calibrated  for  the  desired  positive 
or  negative  range  measuring  a  pressure  of  the  opposite  sign 
will  cause  the  zero  point  of  the  transducer  to  shift.  Re¬ 
calibration  will  then  be  required.  Because  both  positive  and 
negative  differential  pressures  were  measured  it  was  found 
that  one  transducer  had  to  be  calibrated  for  pressures  less 
than  atmospheric  and  one  for  pressures  greater  than  atmos¬ 
pheric.  A  Sample  calibration  program  is  included  in  Table 
C-I. 

Once  the  calibration  process  is  complete  the  data 
acquisition  program  is  run.  This  program  will  record  the 
required  pressures  and  temperatures  and  upon  completion 
store  the  values  under  the  desired  filename  on  the  flexible 
disc.  The  data  reduction  program  is  then  loaded  from  the 
disc  storage  and  executed.  Upon  completion  of  the  data 
reduction  program,  plot  programs  are  called  which  will 
provide  output  on  the  722S-A  Plotter. 
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The  data  acquisition  system  as  presently  configured  is 
semi-automatic.  Temperatures  can  be  automatically  taken  in 
sequence  but  the  program  must  stop  and  allow  the  operator  to 
manually  step  the  scanivalve  and  arrange  the  valving  for  the 
correct  transducer.  The  use  of  a  digital  stepping  device 
would  further  automate  the  process.  Additionally,  a  trans¬ 
ducer  dedicated  to  only  determining  positive  or  negative 
differential  pressure  and  electrically  operated  manifold 
valves  would  complete  the  process. 

Even  though  the  system  is  not  yet  completely  automatic, 
the  ability  to  reduce  data  on  site  provides  more  flexibility 
and  ease  of  conducting  the  thesis  research. 
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430  MAT  PRINT  X 
440  END 


10  1  The  name  of  this  pro?ratt 

is  calssxo 

20  OPTION  BASE  0 

30  OISP  "Ttiis  program  is  desian 
ed  to" 

48  DISP  "calibrate  the  0-35  in 
t  rans-" 

50  DISP  “ducer.lt  will  sive  you 
values" 

60  DISP  "to  input  into  the  CHEB 
V  pro-" 

70  OISP  "a ram  to  set  coefficien 
t  s  tor  a " 

80  OISP  "calibration  function.il 
ake  sure" 

90  DISP  "that  the  55  in  transd 
ucer  is" 

100  DISP  "DISCONNECTED  from  the 
manifold. " 

110  DIM  X < 4 , 1  ) 

120  DISP  "With  the  35  in  transdu 
cer  dis-“ 

130  DISP  "connected  at  the  mani f 
old j  set “ 

148  DISP  "the  zero  point  on  the 
carrier" 

150  DISP  "demodulator . " 

160  OUTPUT  709  s " UR5ARAC3SA" 

170  DISP  "Press  CQHT  when  ready. 

ft 

180  PAUSE 

190  DISP  "Reconnect  transducer  t 
o  the" 

200  DISP  "mani fold . When  air  Pres 
sure  is" 

210  DISP  "set  at  20  in  H20  press 
CONT. “ 

220  PAUSE 
230  1=2© 

240  FOR  J=4  TO  0  STEP  -1 
250  X<J,0>=I 

260  OUTPUT  709  i"0R5AI3SA" 

270  ENTER  709  j  X<J,l> 

280  OUTPUT  709  ;  ''ARAC3" 

290  DISA  X<J,0>,X<J- 1) 

300  1=1-5 

310  OISP  "Decrement  IN  H20  by  5 
IN  to  " ; I 

320  IF  J=0  THEN  340 
330  PAUSE 
348  NEXT  J 

358  DISP  "  P  V" 

360  DISP  "inH20  voltase" 

370  MAT  DISP  X; 

380  DISP  "HARD  COPY<l=yes  2=no>" 

390  INPUT  N 

400  IF  N#1  THEN  440 

410  PRINT  “  P  V" 

420  PRINT  " inH20  voltaee" 


TABLE  C-I.  SAMPLE  CALIBRATION  PROGRAM 
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FIGURE  C-I.  SCHEMATIC  OF  DATA  ACQUISITION  SYSTEM 


FIGURE  C-IIc.  ILLUSTRATION  OF  DATA  ACQUISITION  SYSTEM 
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FIGURE  C-IIb.  ILLUSTRATION  OF  DATA  ACQUISITION  SYSTEM 


APPENDIX:  D 


UNCERTAINTY  ANALYSIS 

The  determination  of  the  uncertainties  in  the  experimen¬ 
tally  determined  pressure  coefficients,  pumping  coefficients 
and  velocity  profiles  was  made  using  the  method  described  by 
Kline  and  McClintock  [6].  The  uncertainties  obtained  by 
Ellin  [1]  using  the  second  order  equation  suggested  by 
Kline  and  McClintock  [6]  are  all  applicable  to  the  experi¬ 
mental  work  reported  herein  and  are  summarized  in  the  follow 
ing  table. 


TABLE  XIV 


UNCERTAINTY  IN  MEASURED  VALUES 


a 

AP 

PV 

P 

u 

APS(+) 
AP  (**) 


or 

AP 


or 


±  1  R 
±  1  R 

±  0.01  psia 
i  0.01  in.  H20 
±  0.01  in.  Hz0 
±  0.05  in.  H20 
±  0.01  in.  H20 
±  0.01  in.  H20 
±  0.01  in.  H20 
±  0.20  in.  H20 
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T 


or 


T 


a 


±  1  R 
±  1  R 


PT  (+++) 


±  0.1  in.  H20 


UNCERTAINTY  IN  CALCULATED  VALUES 


AP* 

T*~ 

W*T* 


1.9% 


1 . 4 1 


V/V 


avg 


(  +  ) 


(++) 


(+++) 


2.5% 

The  pressure  differential  across  the  secondary 
flow  nozzles,  P  ,  is  the  major  source  of  un¬ 
certainty  in  the  pumping  coefficient. 

The  pressure  differential  across  the  tertiary 
flow  nozzles,  P^,  is  the  major  source  of  un¬ 
certainty  in  the  pumping  coefficient. 

The  measurement  of  the  total  pressure  for  the 
velocity  profile  is  the  major  source  of  un¬ 
certainty  in  the  velocity  calculation. 


165 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  0142  2 

Naval  Postgraduate  School 

Monterey,  California  93940 

3.  Department  Chairman,  Code  69  2 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93940 

4.  Professor  Paul  F.  Pucci,  Code  69Pc  5 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93940 

5.  LT  Charles  R.  Ellin  1 

13512  Westwind  Drive 

Silver  Spring,  Maryland  20904 

6.  Mr.  Charles  Miller  1 

NAVSEA  Code  0331 

Naval  Ship  Systems  Command 
Washington,  D.C.  20362 

7.  Mr.  Olin  M.  Pearcy  1 

NSRDC  Code  2833 

Naval  Ship  Research  and  Development  Center 
Annapolis,  Maryland  21402 

8.  Mr.  Mark  Goldberg  1 

NSRDC  Code  2833 

Naval  Ship  Research  and  Development  Center 
Annapolis,  Maryland  21402 

9.  Mr.  Eugene  P.  Wienert  1 


Head,  Combined  Power  and  Gas  Turbine  Branch 
Naval  Ship  Engineering  Center 
Philadelphia,  Pennsylvania  19112 


166 


I  10. 

Mr.  Donald  N.  McCallum 

NAVSEC  Code  6136 

Naval  Ship  Engineering  Center 

Washington,  D.C.  21362 

1 

LT  J.P.  Harrell,  Jr.,  USN 

2004  Cloverleaf  Place 

Ardmore,  Oklahoma  73401 

1 

12. 

LT  J. A.  Hill 

Route  2,  Box  61 

Wood  River,  Nebraska  68883 

1 

13. 

1 

LT  C.M.  Moss 

625  Midway  Road 

Powder  Springs,  Georgia  30073 

1 

14. 

LT  R.J.  Lemke 

2902  No.  Cheyenne 

Tacoma,  Washington  98407 

1 

15. 

LCDR  D.L.  Ryan  III 

6393  Caminito  Luisito 

San  Diego,  California  92111 

2 

16. 

LT.  R.S.  Shaw 

147  Wampee  Curve 

Summerville,  South  Carolina  29483 

1 

' 

17. 

LT.  C.P.  Staehli 

Route  2,  Box  64B 

Burton,  Washington  98013 

1 

18. 

t  ! 

M 

1 

LT.  D.  Welch 

1036  Brestwick  Commons 

Virginia  Beach,  Virginia 

1 

i  19. 

1 

Dean  of  Research,  Code  012 

Naval  Postgraduate  School 

Monterey,  California  93940 

1 

! 

J 

i 

J 

i 

! 

* 

i 

M 

[ 

167 

*91 

END 


DATE 

FILMED 


DTIC 


